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Abstract 

The Quasiparticle Random Phase Approximation equations are 
solved taking into account the Pauli Principle at the expectation value 
level, and allowing changes in the mean field occupation numbers to 
minimize the energy while having the correct number of particles in 
the correlated vacuum. The study of Fermi pn excitations in 76 Ge us- 
ing a realistic Hilbert space shows that the pairing energy gaps in the 
modified mean field are diminished up to one half of the experimental 
value when strong proton-neutron correlations are present. Addition- 
ally, the Ikeda sum rule for Fermi transitions is violated due to the 
lack of scattering terms in the phonon operators. These results call for 
a critical revision of the double beta decay half-lives estimated using 
the QRPA extensions when standard QRPA calculations collapse. 
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1 Introduction 



The random phase approximation (RPA) and its quasiparticle generalization 
(QRPA) have been widely used in the last decades to study electromag- 
netic transitions and beta decays in medium and heavy nuclei Jl], |2|]. The 
proton-neutron quasiparticle random phase approximation (pn-QRPA) has 
been extensively employed in the description of single and double beta de- 
cays in vibrational nuclei. However the RPA develops a collapse, i.e. it 
presents imaginary eigenvalues for strengths beyond a critical value of the 
force § 1 1, |. 

A whole family of extensions of the QRPA, called renormalized QRPA 
(RQRPA) are known that do not develop any collapse by implementing the 
Pauli principle in a consistent way, beyond the simplest quasiboson approx- 
imation |7|, |8|, |9|, I0[ 11 , 12, [13|] . However, in its simplest versions there is a 



violation of the non energy weighted Ikeda sum rule 0]. Calculations to de- 
termine the amount of the violation and some improvements to the RQRPA, 



in order to restore the sum rule, have been presented fl5fl . It has been shown 
that treating simultaneously BCS and QRPA equations one can fulfil the 
Ikeda sum rule for the Fermi case when a schematic model is used |16| . 



In recent articles we studied the expectation values of the particle and 
quasiparticle numbers, the particle number fluctuations and the number of 
particle pairs with J=0, T=l and J=l, T=0 in the ground state of 76 Ge as 
a function of the residual proton-neutron interaction, using realistic Hilbert 
spaces [|T|], |18| . We found an important amount of particle number fluctua- 
tions in the RQRPA ground state beyond the QRPA collapse, pointing out 
a source of uncertainty in the RQRPA results. The analysis of the number 
of pairs showed that the isoscalar-isovector phase transition found in exact 
calculations is causing the QRPA collapse and is missed in the RQRPA for- 
malism. 

In the present work we go a step further, studying 76 Ge with a renor- 
malized RPA where at the same time the mean field is changed by minimiz- 
ing the energy and fixing the number of particles in the correlated ground 
state. While particle number fluctuations are smaller than in the previous 
cases, they still exhibit a clear increase after the point of collapse. More 
remarkably, the pairing gap is strongly reduced in comparison with its ex- 
perimental value, and the Ikeda sum rule is violated. Both results cast serious 
doubts about the double beta decay half-lives estimated using the QRPA ex- 
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tensions, in particular for those nuclei where standard QRPA calculations 
collapse H, |l(], [12| . Together with the conclusions obtained in |17|, 0, we 
get a clear picture of the limitations associated with the QRPA extensions, 
which at the end are the same common sense ask for: you cannot allow the 
residual proton-neutron interaction to dictate the composition of the ground 
state wave function without missing contact with actual nuclei, even if the 
formalism allows you to overcome the collapse. 

In Section 2 the renormalized gap and number equations are introduced, 
whose relationship with the RQRPA equations is shown in Section 3. Some 
relevant expectation values and the Ikeda sum rule are discussed in Section 
4, results for 76 Ge are presented in Section 5 and the conclusions in Section 
6. 



2 Renormalized Gap and Number Equations 

In this section the gap and number equations are obtained minimizing the 
Hamiltonian expectation value (0|i?|0) in the correlated vacuum |0). The 
Hamiltonian is 

H = H p + H n + H pn . (1) 
The first two terms refer to the proton and neutron hamiltonians 

H t = Y,( e t ~ ^) a t a t + \ < t i t 2|^|t 3 t 4 > a^a^e^ata, (2) 

t 1 t's 

where the single particle energies are denoted by e t , the chemical potential 
by A, and the last term corresponds to the proton-neutron interaction 

H P ,n= < p,n|V|p',n' > aJ,aj/v<V • ( 3 ) 

p,p',n,n' 

The subscripts t(t) stand for p(p) (protons) or n(n) (neutrons), being t = 
t, m t , with t = {n t , l t ,jt} and m t = rrij. 
Through the Bogoliubov transformation 

"t = u t<A - v t ai, (4) 
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with aj = (-l) t+mt a t - mU we get [|TJ 
H = U + J2 Ztopepftp + E ™ n ejf n + H 22 + H m + H 04 + H 13 + H 3U (5) 



being 



Nut = 



(2jt + 1) 



and where C/ and the quasiparticle energies et are defined as 



^ = E 



2Q t v?(e t - A) + E Jn t ,n t vtvlF(tt, t't', 0) - JWtA, 



(6) 



being 



t> V "t 



'a, 



(v 2 - M?) + 2A t U t V t 



At = -1/2 E Jtt u v v v G{tt,t't',G), 



(7) 



(8) 



the 'gap' and F, G the usual particle-hole (PH) and particle-particle (PP) cou- 
pled two-particle matrix elements. The terms H nm in (H) destroy m quasi- 
particles and creates n quasiparticles, respectively. 

To obtain the quasiparticle mean field occupations t>j the ground state 
energy (0|if|0) is minimized, i.e. it is asked that Q 



d_ 

dv 



(0\H\0) = 0, 



subject to the constrictions 

(0|JV|0) = iV, <0|Z|0) = Z, 



(9) 



(10) 



and the normalization uf + vf = 1, being N and Z the neutron and proton 
particle number operators respectively. 
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The standard BCS procedure treats protons and neutrons separately and 
in Hamiltonian (||) only pairing interactions, refering in general to like par- 
ticles interacting through the J=0 channel, are included . In this case the 
ground state energy is just U as defined in Eq. @ and |0) is the BCS ground 
state. The residual interactions, either between like particles or protons and 
neutrons, represented by the terms H nm in (|5]), are usually included in a 
second step, most commonly using the QRPA [|l|, [^]. 

One of the drawbacks of this procedure is that Eq. (|T0|) is only enforced 
for the BCS vacuum. When residual interactions are present the expecta- 
tion values of the particle number operators do not coincide with the actual 



number. The SCRPA pOj , [11] , f21|j is designed to overcome this difficulty by 
solving Eq (Q) and ([TO]) using the RPA vacuum. 

Full selfconsistency requires to consider the proton-neutron interaction 
contribution in the minimization, but the system of equations which de- 
scribes this problem is nonlinear and rather complicated, and has only be 



implemented in schematic models |L1, |22|, [23] . In order to perform calcu- 



lations in realistic Hilbert spaces we will include in this first step only the 
like-particles part of the Hamiltonian, while the ground state |0) will be sen- 
sitive to proton-neutron interaction through the RPA equations, following a 
philosophy close to Ref J24|. From here on we will refer to this approxima- 
tion as SRQRPA. We are absolutely aware that our treatment is not fully 
self-consistent. We are just meeting the requirements of Eq. ([10|) for the 
RPA vacuum, and including the modifications in the mean field due to pn- 
correlations at the lowest level. However, as is shown below, even this mild 
modifications have very important effects in the observables of the system. 



Following [pi[| , but including only the modifications of the gap equations 
due to use of the RPA vacuum instead of the BCS vacuum, i.e. not taking 
into account the proton-neutron residual interaction explicitly, we arrive to 
the modified gap equation 



2(e t - \)v t u t 



[u. 



(11) 



where 



e t = e t + J2\l-^-4F(tt,t't',0)) 



(12) 
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is the single particle energy corrected by the self-energy and 



A t = -l/2j2\^u tl v t ,G(tt,t , t',0)(l-2(0\Af t ,\0)), (13) 

t , V "t 

the 'renormalized' gap. Notice that the pairing interaction appears in the 
gap equation renormalized by the presence of proton-neutron residual inter- 
actions, which introduce a finite number of quasiparticles in the RPA vacuum 
[|T7|]. For the sake of simplicity we have dropped some higher order terms in 
Eq. (|TTD which come from the interaction between like particles connect- 
ing the many quasiparticle components of the correlated ground state. The 



definition of the renormalized gap ([13]) is unique at this level of approxima- 



tion This is an important definition to be kept in mind, because the 
factor (1 — 2(0|JV t '|0}) will play a definite role in suppressing the gap when 
proton-neutron residual interactions are large. 

What remains is to couple the renormalized gap problem with the RQRPA 
equations. 



3 RQRPA 

The nuclear excited states are constructed as || |15 



\\JM) = n\\JM) |0), (14) 
fit(AJM) = Y, [x pn (XJ)Al n (JM) - Y pn (XJ)A pn (7M)] , (15) 



pn 

where 



AUJM) = [ala&) JM D£' 2 , D pn = (l - (0\Af p + Af n \0)) , (16) 



are the renormalized two-quasiparticle proton-neutron creation operators, 
which satisfy 

(0| [A pn (JM), Al, nl {J'M')\ |0> = 5 pp ,5 nn ,5jj,5 MM , (17) 
Here |0) is the RPA correlated ground state, defined by the condition 

n(XJM)\0}=0. (18) 
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Each amplitudes X pn (XJ) and Y pn (XJ) are associated with the excitation 
energy uo\j of the A-th state with angular momentum J. They are the eigen- 
vectors and eigenvalues, respectively, of the RPA equations |TJ, |2j 

( A(J) B(J) \(X(XJ)\_ ( X(XJ) \ 

[ B*(J) A* {J) J { Y(XJ) J XJ { -Y(XJ) J ' 1 ' 



where 



B pn ,p'n'{J) = Dl^Ul L p / n /(J)-DpC> (20) 



with 



U pn, P 'n'( J ) = G(pn, pn,JM) {u p u n u p >u n > + v p v n v p ,v r , 
+ F(pn, p'n', JM) (u p v n u p >v n i + v p u n v p >u r , 



U pn, P 'n'( J ) = -G(pn,p'n',JM)(u p u n v p ,v n , + v p v n u p >u n ,) 

+ F(pn,p'n',JM)(v p u n u p iv n >+u p v n v p >u n >). (21) 

The RPA equations ( jT3 ) depend explicitly on the mean field occupations 
Vi through the A and B matrices defined in (|20|) and (21). At the same 



time, the quasiparticle occupations Vi depend on the RPA amplitudes X, Y, 
on which depend also |0), through the renormalized gap equation ( ]TT|) and 
number equations fllPl) . 

4 Expectation values in the correlated RPA 
vacuum 

Using the quasi-boson approximation 



A pn (JM),Al, n ,(J'M')\ » (0| [A pn (JM),Al, n ,(J'M')\ |0) = V<W£/J'<W : 

(22) 

the RPA ground state defined by Eq.(^) can be written as 

|0> = N e S \BCS), (23) 
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with the BCS vacuum defined by the property 

a t \BCS) = 0, (24) 

and 

S=l £ _ yj (2 J + 1)C( J) pnp , n , [Al( J)A\, nl { J)] ° , (25) 

C(J) pnp > n > = ^ Y ( J) ^ X X ( J) * x y n , . 



2 

pnp'n'J 



The quasiparticle occupations (0|jVt|0) evaluated using the QRPA vac- 
uum (E3I) have the explicit form 



XJn' Zi L P 

(0|^n|0) = Y, { -^^- D P>n\Y(J) P >n,x\ 2 , (26) 
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The mean particle numbers of protons and neutrons are |25 

(Z) = (0\Z\0) = 2Y,^ P v 2 p + 2^ P (ul-v 2 p )(0\Af p \0), (27) 

p p 

(N) = (0\Af\0) = 2 £ ^ + 2 £ - t£) (0|Ago>, 

n n 

When the BCS vacuum is used, the second term in each of these expressions 
vanish and one is left with the usual number equation. When the residual 
interaction is present these terms have a relevant contribution to the particle 
number JT7|] . They are included in the present work in a selfconsistent way 
by modifying the mean field occupations v p , v n to obtain the correct number 
of particles in the correlated RPA vacuum. 

The Fermi transition operators are written in terms of the pair creation 
and annihilation operators as [§. [| 



r ~ = EK tt «[«!«if° + u p v n [a p a n ] 00 ] , r + = {r'}^ (28) 

pn 
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which are not the exact expressions because they are missing the scattering 
terms which create a proton (neutron) quasiparticle and annihilates a neutron 
(proton) quasiparticle. 

The total strengths S± associated with these transition operators are 

5 ± = 53l<AJ=0|T i |0>| a . (29) 

A 

The Ikeda sum rule states that, when the exact operators are used, and 
when the set of states |A) is a complete one, including all the states in the 
odd-odd nuclei which can be connected with the ground state |0) of the 
even-even nuclei through the transition operators, then 

S_ - S+ = N - Z (30) 

However, in the present case the Fermi transition operators are truncated, 
and the strengths difference has a more complicated form |25[] 



S_ - S + = (0\N - Z\0) - - v 2 p )(0\Af n - Ago). (31) 



pn 



The first term is equivalent to N — Z due to the constrictions (|10|). The 
second term gives rise to the violation of the Ikeda sum rule (|3T)D. When the 
BCS vacuum is employed, these term has no contribution. For this reason the 
standard pn-QRPA fulfils the Ikeda sum rule. There are some special cases 
in which the proton and neutron quasiparticle occupations in the ground 
state are equal, and this term also vanishes, as it was found in the SO (5) 
model |16| and in the single mode model of the RQRPA |T5[ . Only when 
the sum in (|^) is restricted to J = 0, forcing j p = j n and Q p = Q n , it 
is true that (0\Af n — Ap|0) = for each level. In the present case, when a 
realistic Hilbert space is employed and the sum over different J's is included 
in (^), those occupations are not equal and the Ikeda sum rule is violated. 
It is worth to mention that, while in general the quasiparticle occupations 
for protons and neutrons are different, the expectation value of the total 
number of quasiparticles (N qp ) = St(0|-^ v t|0) is the same, as can be easily 
seen making the sum in Eq. (I 
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5 Results 



In the present work we study Fermi beta excitations in 76 Ge. We adopt a 8- 
type residual interaction already used previously fL7|, an d our Hilbert space 



has six single particle energy levels, including all the single-particle orbitals 
from oscillator shells 3hu plus lgg/2 and lg?/2 from the 4hu oscillator shell. 
They were obtained using a Coulomb-corrected Wood-Saxon potential. Their 
numerical values for 76 Ge are tabulated in Table 1 of ref. |26|. We include 
J 71 = ± , .., 3 ± in the sums in Eq.( 



In order to describe the dependence of the various observables on the 
proton-neutron residual interaction we use the parameter 

v pp 
s = pn 

ypair ' 

which is the ratio between the coupling constant in the proton-neutron J = 
particle-particle channel and the pairing force constant v pair = (fpp ir + 
v P ^ r )/2. We compare results obtained within the usual QRPA with those 
coming from the RQRPA (no mean field modification), and the SRQRPA. 

In Fig. 1 we show the lowest J n = + excitation energies u (in MeV) in 
76 As, calculated from the 76 Ge ground state, for the QRPA (dotted lines), 
RQRPA (dashed lines) , and SRQRPA (full lines) approximations, as a func- 
tion of the residual interaction parameter s. It can be clearly seen that the 
excitation energy goes to zero, and collapses around s = 2.5 for the QRPA, 
while in the RQRPA and SRQRPA formalisms the collapse is avoided, al- 
though the excitation energies predicted after the collapse in these two ap- 
proaches differ by a factor two. 

In Figure 2 the sum rule — S + , normalized to the Ikeda value N — Z, 
is presented as a function of s with the same convention of Fig. 1. While 
in the QRPA the Ikeda sum rule is always fulfiled, in the other cases it is 
violated. Notice that, having the correct number of protons and neutrons in 
the correlated ground state, the SRQRPA departs less from N — Z that the 
RQRPA, but the departure is anyway noticeably. It is important to point out 
that, in the present calculations, even when s = the other proton- neutron 
residual interaction channels, for J 7^ 0, are present and have a finite value. 
For this reason there is a violation of the Ikeda sum rule at s = 0. 

To describe in more detail the above mentioned point we present in Figure 
3 the total number of quasiparticles as a function of s. When only J = 
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states are included in the sum in Eq. fl26|) the number of quasiparticles goes 
to zero when s is zero. In the other three cases the residual interaction in 
the other channels is present and increase this number. 



The fluctuations in the particle number AN = y (0|(iV — iV) 2 |0) repre- 
sent a warning about the region of applicability of the different QRPA models 



17], [18j . They are shown in Fig. 4a (for protons) and Fig. 4b (for neutrons) 



as a function of s. It is interesting to note that the pure BCS fluctuations 



AN = 2yY^t ^t u t v t diminish when the pn residual interaction increases, due 
to the sharpening in the particle distribution when the mean field is varied 
to keep the number of particles in their right value. In the SRQRPA these 
fluctuations diminish slightly at first, dominated by this mean field effect, to 
increase in a noticeably way after the QRPA collapse. The QRPA results for 
the total quasiparticle number and the particle number fluctuations explode 
for s pa 2.5 when only the J n = + contribution is included in Eq. (|26|). 
When all the angular momentum contributions are taken into account the 
J w = 1 + component collapses around s = 2.3, and these expectation values 
diverge. 

The single particle occupation numbers Vp,v^ are presented in Figure 5 
as a function of the single particle energies. While only six levels are in use, 
the curves represent their smoothed values. It can be seen that both protons 
and neutrons have their distributions sharpened, i.e. the larger s the more 
they resemble a step distribution, although this effect is more evident for 
protons. 

The behavior of the pairing energy gap as a function of s is shown in Fig. 
6a (for protons) and Fig. 6b (for neutrons). As advanced in the previous 
section, the pairing gaps, while different for the different single particle levels, 
all behave the same way. They show a strong reduction when s increases, 
down to one half of the original value. This is not a curiosity: the pairing 
energy gap is observable, and any attempt to properly describe the beta 
decay using the QRPA extensions must keep the descriptive power of the 
simplest formalism, the QRPA, where the pairing constants are adjusted to 
reproduce the observed gap. 
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6 Conclusions 



A study of the beta decay of 76 Ge using the SRQRPA was presented. The 
mean field gap and number equations were introduced and solved together 
with the RQRPA equations, but the selfconsistency was only included to 
guarantee the good number of particles in the correlated ground state. 

The Ikeda sum rule was studied in some detail. It is known to be fulfiled 
in standard pn-QRPA calculations ||, |]], violated in RQRPA ones [14| and 
recovered when self-consistent RQRPA calculations are performed in simple 



models [|T5|, |16 |. In the present work it was shown that the Ikeda sum rule is 



violated when a realistic Hilbert space is used in spite of using the number 
selfconsistent QRPA approach. 

The most remarkable phenomenon found in the calculation was the strong 
reduction of the pairing gap. In the presence of the proton- neutron residual 
interaction the mean field changes needed to have the correct number of 
particles generate sharper distributions for the occupations of single particle 
states, visible when plotting v 2 as a function of the single particle energies 
for different values of s. At the same time, it strongly reduces the energy 
gap, which can be as small as a half of the observed value. 

Having in mind the various merits of the renormalized and self-consistent 
extensions of the quasiparticle random phase approximation, we conclude 
that it must be used with extreme caution in the region where the stan- 
dard QRPA collapse. Not only the particle number fluctuations have a clear 
increase, but observable quantities depart from their measured values. 
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Figures caption 



Figure 1. Lowest J w = + excitation energies in 76 As, calculated from the 
76 Ge ground state, for the QRPA, RQRPA, and SRQRPA approximations, 
as function of the residual interaction parameter s. 

Figure 2. Fermi sum rule S- — S+, normalized to the Ikeda value N — Z, 
as function of s with the same line conventions used in fig. 1. 

Figure 3. Total quasiparticle number (equal for protons and neutrons) 
evolution as function of s. The result for the QRPA case when the sums in 
Eq. (26) are restricted to J 71 = + is also shown. 

Figure 4. Proton (a) and neutron (b) number fluctuations obtained using 
the BCS, QRPA with only J n = 0+ and all the J's in Eq. ([§), RQRPA and 
SRQRPA descriptions. 

Figure 5. Single particle occupation numbers v% (thick lines), and v\ (thin 
lines) as a smoothed function of the single particle energies e, for the values 
s = 0., 2.5, and 3.1 of the residual interaction parameter. 

Figure 6. Evolution of the pairing energy gap as a function of s for protons 
(a) and for neutrons (b). 
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